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(57) The present invention comprises a simulator 
including at least one data storage apparatus 
(71) and at least one processor (72). The 
simulator further has a material property data 
group (44), a system state data group (42) com- 
prising geometrical data for defining a plurality 
of geometrical subvolumes and system condi- 
tion data representing the state of each sub- 
volume of a simulated system, an event 
probability data group (46) including a plurality 
of event subgroups for uniquely identifying a 
set of events each having a probability of oc- 
currence represented by the event probability, 
each of the event probabilities being associated 
with at least one of the geometrical sub- 
volumes, and an event process data group (48) 
including the time dependent process functions 
for all events which can occur in the system. 
These data groups and at least one functional- 
relationship link for correlating all of them, 
whereby the material energy system can be 
represented by such data groups without re- 
quiring an explicit multiple dimension data 
structure, are stored in the data storage ap- 
paratus. The processor includes a data input 
means (73) to read these data groups from the 
data storage apparatus. The processor further 
has means for performing a simulating process 
comprising the steps of: (a) computing the 
probability of occurrence for each of the events 
and determining a time step by using a total 
probability of occurrence as a sum of the prob- 
ability of occurrence for each of the events ; (b) 
randomly selecting a simulated event for the 
time step and updating the system condition 
data group representing the system state result- 
ing from the occurrence of the selected simu- 
lated event ; and (c) repeating the steps (a) and 
(b) until a predefined simulation end time is 
reached. 
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This invention relates generally to a simulator for 
simulating, predicting, and/or controlling a time- 
progressive material-energy system. More particular- 
ly, this invention teaches a stochastic simulator to 
perform real-time simulation, prediction or control of 5 
a spatially inhomogeneous system by fully utilizing 
the data storage capacity and processing power pro- 
vided by an on-line computer system. 

The usefulness of a simulator for a highly dynam- 
ic and complex system is often limited to design or 10 
training but seldom used for real time control. There 
is an intrinsic limitation, namely, the processing time 
required to propagate the system by a simulator in 
many cases is slower than the rate of occurrence of 
the real events taking place in the simulated system is 
thus preventing the simulators from being used on a 
real time basis. 

The usefulness of a conventional simulator is also 
limited by the fact that a simulator in most cases 
is 'system specific', i.e. a simulator is constructed to 20 
only simulate a system with certain geometrical con- 
figurations, material compositions, physical or chem- 
ical processes, and dynamic characteristics with very 
limited ranges of flexibility. A simulator originally de- 
signed to be implemented on a certain combination of 25 
hardware and software for simulating a specific sys- 
tem is generally very rigid in its data structures and 
processes of operations so that the simulator can not 
easily be converted to either simulate a different sys- 
tem or be implemented on a different combination of 30 
hardware and software architectures. With the wide- 
spread use of the more intelligent processors and the 
increased capacity of data storage and computer net- 
works to enhance the speed and performance of dis- 
tributed parallel and concurrent data processing, rigid 35 
data structures and inflexible simulation processes 
pose a major hindrance preventing the simulator art 
from taking advantage of the ever increased data 
processing power and data storage capacity provided 
by the modern electronic and computer industry. 40 

Simulating and predicting the state of a time- 
dependent material-energy system can be a very dif- 
ficult and time consuming task for several reasons. 
Traditionally, a mathematical formulation which often 
involves a set of complicated and dynamically interac- 45 
tive differential equations is used to explicitly repre- 
sent various space and time-dependent processes 
which may occur under a variety of system conditions. 
Explicitly means that all the information predeter- 
mined to be necessary to represent the entire system so 
is actively determined for each time step for the sim- 
ulation. 

To do this, the simulator performs a process to 
solve this set of differential equations by using a finite 
difference or finite element method to obtain the nu- 55 
merical solutions. For complex systems in particular, 
simulation and prediction inaccuracies may result 
from either the approximation errors caused by the 



finite difference or finite element solution methodol- 
ogies or from the fact the differential equations may 
not accurately represent the system behavior at all 
times in all locations over the entire system. 

Instead of treating the simulation of a time- 
dependent material-energy system as a mathemati- 
cal solution process, stochastic simulators have been 
developed. These operate by continuously and dy- 
namically tracing a set of individual events which may 
occur over the time-span of simulation. The simulator 
makes an assumption that each event has a certain 
probability of occurrence which is determined by the 
system conditions and material properties at a specif- 
ic point of time. In a predefined time-window which is 
commonly referred to as a time step, any one of the 
simulated events may occur. There is, on the aver- 
age, one occurrence of an event within a time step 
when the length of that time step is selected to be 
equivalent to the inverse of the total probability of oc- 
currence of all events. The likelihood of occurrence of 
certain event within that time step is proportional to 
the probability of occurrence of that individual event. 
A stochastic simulator commonly comprises at least 
one processor which is capable of defining a time step 
and then randomly selecting an event as a triggering- 
event for that time step. The likelihood of an event be- 
ing randomly selected is proportional to the probabil- 
ity of occurrence of the selected event relative to that 
of the other events. As the consequence of the occur- 
rence of the selected event, the processor updates 
the system conditions by use of material property ta- 
bles, system process equations and other required 
data groups stored in at least one data storage appa- 
ratus. The processor then selects a new time step and 
repeats the event selection and system condition up- 
date processes to propagate a system through the 
entire time span for system simulation. 

A stochastic simulator was first used to simulate 
a spatially homogeneous chemical system where 
many chemical reactions occur (D. L. Bunker et al., 
Combustion and Flame, Vol. 23, p.373 ff., 1974). The 
probability of an individual chemical reaction is de- 
pendent on the concentrations of a variety of chemi- 
cal reactants. The functional dependence and the 
concentrations are stored in the storage apparatus. 
The processor generates a random number to select 
among various chemical reactions a specific reaction 
in each time step thus allowing the simulation to prop- 
agate in time. 

The time base for determining the length of time 
step is based on the total sum of all reaction probabil- 
ities. Two specific embodiments of this type of simu- 
lator are described in documents available from 
Quantum Chemistry Program Exchange (Indiana Uni- 
versity) as Simulator 293 and Simulator QCMP069. 

The tasks performed by a simulator according to 
Bunker et al. are divided into many stages as illustrat- 
ed in Appendix A where a stochastic simulator is used 
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to simulate a spatially homogeneous chemical reactor 
10 The chemical reactor 10 may be described as a 
very small volume containing a limited but adequate 
number of molecules 12 of the various species in- 
volved in a set of chemical reactions. The initial state 
of this volume, e.g. temperature, pressure and con- 
centrations, is stored in the data storage apparatus of 
the simulator (see Appendix A, step (b) of stage 1). 
The intrinsic characteristics of the system, e.g. the al- 
lowed reaction pathways and their rate of occurrence 
(step (a)) , are also stored in the storage apparatus of 
the simulator as a plurality of system process arrays. 
The processor first computes the initial probabilities 
from the initial state and the system process arrays. 
The processor then propagates the simulation by sto- 
chastically selecting one event at a time and iterative- 
ly executing the sequence of steps listed in stage 3. 
The simulator then generates a data group for storage 
in the data storage apparatus in stage 4 representing 
the state of the chemical reactor at each time step 
when the processor completes the simulation steps 
listed in stage 3. 

Figs. 2A, 2B and 2C schematically illustrate the 
event selection operation used by Bunker et al. for 
simulating the chemical reactor 10 which comprises 
three chemical reactions, i.e. reactions 1, 2 and 3, 
each has a reaction rate of R[1], R[2] and R[3] respec- 
tively, defined by rate constants k[1], k[2] and k[3] and 
the instantaneous concentrations of the species A, B, 
and C. The rates R[n] can be considered to be the 
probabilities of reactions 1 to 3 occurring per second, 
and the total probability P (per second) is then the 
sum R[1] + R[2] + R[3]. Since the mean probability of 
an event occurring in the time period 1/P is unity, that 
period is taken as the time step for a single event and 
the elapsed time in the simulation is advanced by that 
amount (Fig. 2B). The processor uses a random num- 
ber to select one of the three reactions. Depending on 
the value of the random number any of the three re- 
actions can be selected (Fig. 2C), with the likelihood 
for selection of a given reaction equal to its contribu- 
tion to the total probability P. This is represented by 
the lengths of the line segments in Fig. 2C. 

Another chemical simulator has been developed 
by Gillespie wherein the simulator generates two ran- 
dom numbers; one for identifying a chemical reaction 
and another for properly weighting the time base 
used for time step determination. (Journal of Compu- 
tational Physics, Vol. 22, No. 4, pp. 403-434, 1976). 
Better accuracy is achieved because the time step is 
dynamically weighted by the probability of occur- 
rence for each individual event. 

The stochastic simulators disclosed by Bunker 
and Gillespie are both homogeneous simulators 
which can only be used to simulate spatially homoge- 
neous systems. Gillespie proposes to generalize the 
application of these stochastic simulators to spatially 
inhomogeneous chemical reactions by first dividing 



the simulated system into a plurality of subvolumes 
and allowing the chemical reactions to proceed "inde- 
pendently" in several subvolumes. The mass trans- 

5 port of reactants and reaction products between the 
subvolumes are proposed to be assigned with prob- 
abilities in a manner similar to that used for the indi- 
vidual chemical reaction steps. A 3-dimensional prob- 
ability function has to be used, reflecting both a reac- 

10 tion step and a physical location. The proposed sim- 
ulator is required to stochastically select a specific 
subvolume as well as a specific reaction step to prop- 
agate the chemical reaction throughout the entire 
system in time. 

15 A simulator as proposed by Gillespie thus would 

be required to simultaneously process three random 
selection processes; one for the selection of subvo- 
lumes, one for reaction type and one for time step 
weighting. Since the proposed simulator selects the 

20 reaction step independently from the selection of a 
particular subvolume, the time step obtained after the 
random time-step weighting process may be con- 
strained by the selected subvolume and therefore 
may not be appropriate for other subvolumes wherein 

25 other reaction steps may have rates of reactions sig- 
nificantly different than that of the selected volume. 
The proposed simulator would thus have difficulties 
in determining a common time base for defining a 
time step for the entire system. Recognizing this lim- 

30 itation, Gillespie states that application of the pro- 
posed simulator to inhomogeneous chemical systems 
is limited to tubular systems with very similar subvo- 
lumes where differences in time evolution of the va- 
rious subvolumes would be small. 

35 Various discrete simulation technologies for 

chemical kinetics are reviewed by J. Turner, Journal 
of Physical Chemistry, Vol. 81, No. 25, pp. 2379- 
2408, 1978). Turner briefly discusses the space- 
dependent stochastic simulation methods with a spe- 

40 erf ic reference to the simulator proposed by Gillespie 
and proposes a method which is essentially an exten- 
sion of a one-dimensional system similar to that of Gil- 
lespie. It is obvious to those of ordinary skill in the art 
that a spatially inhomogeneous stochastic simulator 

45 may have broad applications. However, it is well rec- 
ognized that due to the complexity of the system and 
the need to simultaneously process random selection 
of multiple parameters to accurately trace the 
changes in a dynamic system, a three dimensional 

so spatially inhomogeneous stochastic simulator would 
be too complicated and too slow to be practically use- 
ful. 

It is therefore an object of the present invention 
to provide a stochastic simulation system capable of 
55 accurately simulating, predicting or controlling a simu- 
lated system which is spatially inhomogeneous. 

According to the invention there is provided a 
simulator system for simulating and/or controlling a 
time-progressive material-energy system comprising 
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at least one data storage apparatus including a data 
correlating means; a data structure stored in said 
data storage apparatus divided into data groups de- 
fining geometric configuration, system state, material 5 
properties, material-energy processes, and probabil- 
ities of occurrence of said material- energy processes; 
said data correlating means being capable of defining 
at least one functional-relationship link between data 
stored in said data groups. to 

There is further provided a method of providing 
a data structure for a spatially inhomogeneous sto- 
chastic simulator comprising the steps of generating 
system state data including a plurality of subvolumes 
each identified by a subvolume index (INDX), subvo- 15 
lume system state data and at least one material prop- 
erty identification (ID); generating material property 
data including the properties of the materials con- 
tained in said plurality of subvolumes; generating 
event process data for defining a plurality of events 20 
each having a function for defining a time depend- 
ence of said event as a function of said system state 
for each said subvolume identified by said subvolume 
INDX; generating event probability data for defining a 
probability of occurrence for each of said events in 25 
each said subvolume depending on said system state 
and said material property in said subvolume; and 
storing said system state data, said material property 
data, said event process data and said event probabil- 
ity data in a data storage apparatus wherein said sub- 30 
volume INDXs and said material property identifica- 
tion (ID) provide a correspondence between all of 
said data. 

In a preferred embodiment a plurality of data 
groups stored in said data storage apparatus include 35 
a plurality of data for defining geometric configura- 
tion, system state, material properties, material-ener- 
gy processes and probabilities of occurrence of the 
material-energy processes. These data groups may 
take the form of an array, a linked list, a lookup table 40 
or other types of data structures. The data correlating 
means is capable of defining at least one functional- 
relationship link for correlating all of the data groups 
whereby the material-energy system can be repre- 
sented by the data groups without requiring an explic- 45 
it multiple dimension data structure. The simulator 
further comprises at least one processor including a 
data input means for receiving data of the data groups 
from the data storage apparatus and the processor 
further has means for performing a simulating proc- so 
ess utilizing the functional-relationship link among a 
plurality of said data groups comprising the steps of: 

(a) computing an event probability for each of the pos- 
sible material-energy processes and determining a 
time step by using a total probability of occurrence by 55 
summing up all of said material-energy probabilities; 

(b) randomly selecting a material-energy process 
event for the time step and using the material-energy 
processes, the material property data and the geo- 



metric configuration data for updating the system 
state data resulting from the occurrence of the select- 
ed material-energy process; (c) repeating said steps 
(a) and (b) until a predefined simulation end time is 
reached. The processor further has a control means 
capable of utilizing the updated system state as ob- 
tained from the simulating process in step (b) for con- 
trolling the time-progressive material-energy system. 

In order that the invention may be well under- 
stood a preferred embodiment thereof will now be de- 
scribed with reference to the accompany drawings, in 
which:- 

Fig. 1 is a perspective view of a homogeneous 
material-energy system containing a number of 
molecules of different species; 
Figs. 2A, 2B and 2C describe the event selection 
operation for a homogeneous chemical reaction 
wherein Fig. 2A lists reaction steps and rate laws, 
Fig. 2B lists probabilities and elapsed time calcu- 
lations and Fig. 2C illustrates event selection us- 
ing a random number; 

Figs. 3A and 3B illustrate the data structure in- 
cluding four data groups wherein a subvolume in- 
dex (INDX) and a material property identification 
(ID) provide a link between these data groups; 
Fig. 4 shows a three dimensional material energy 
system comprising many subvolumes; 
Fig. 5 shows a schematic diagram of a soft- 
ware/hardware on-line control system utilizing a 
simulator according to the present invention; 
Fig. 6A lists the probability arrays established in 
the simulator according to the invention; 
Fig. 6B is a schematic illustration of the event se- 
lection operation; 

Fig. 7 shows a three-dimensional material-ener- 
gy system divided into a 3x3x3 array of cells 
wherein a table and conversion equations are in- 
cluded for illustrating the correspondence be- 
tween cell pairs and probability array member, 
said correspondence being an example of how 
INDX may be used to link subvolumes to event 
probability subgroups; 

Fig. 8A shows a chemical reactor controlled by a 
simulator system; 

Fig. 8B shows another chemical reactor control- 
led by a simulator controller to achieve a prede- 
fined reaction rate; 

Fig. 9 shows a material energy system receiving 
energy input from a laser beam; 
Fig. 10 shows the material energy system of Fig. 
9 controlled by a on-line simulator according to 
the present invention; 

Fig. 11 is a comparison of a time-dependent tem- 
perature calculation made with a preferred em- 
bodiment of this invention to a conventional cal- 
culation; 

Fig. 12 is a comparison of a spatial temperature 
distribution calculated with a preferred embodi- 



7 



EP 0 567 302 A2 



8 



ment of this invention to a conventional calcula- 
tion; and 

Fig. 13 is an example of spatial temperature dis- 
tributions calculated for 2-layer systems using a 5 
preferred embodiment of this invention. 
Fig. 3Ashows a data structure 40 comprising four 
data groups. Namely, a system state data group 42, 
a material property data group 44, an event probabil- 
ity group 46 and an event process data group 48. The 10 
system state data group 42 first defines a plurality of 
subvolumes 62 (Fig. 4) each having a subvolume in- 
dex (INDX). These subvolumes represent the geo- 
metrical configuration of a materia I- energy system as 
shown by Fig. 4 where a material-energy system 60 15 
includes a plurality of subvolumes 62. The system 
state data group 42 further provides information of 
the system condition in each of the subvolumes 62 
identified by the associated subvolume index (INDX). 

Each of the subvolumes 62 as defined by the cor- 20 
responding subvolume INDX also contains a plurality 
of materials wherein the properties of each material 
is identified by a material property identification (ID). 
Referring back to Fig. 3, the material property data 
group 44 includes a plurality of data subgroups, each 25 
subgroup defines the material properties identified 
by a material property ID contained in at least one of 
the subvolumes 62 identified by the associated sub- 
volume INDX. 

The event probability data group 46 includes a 30 
plurality of event numbers for uniquely identifying a 
set of events each associated with at least one sub- 
volume 62 identified by the subvolume INDX and 
each having an event probability depending on the 
system condition and the material properties in the 35 
associated subvolumes 62. The event process data 
group 48 includes a plurality of process functions for 
determining the event probabilities from the system 
conditions as provided by the system state data group 
42 in each of the subvolumes identified by an asso- 40 
elated subvolume INDX, and by the material property 
data group 44. 

Specifically, each of these data groups in a pre- 
ferred embodiment is described as follows: 

1. The system state data group 42 first defines a 45 
plurality of subvolumes which represent the geo- 
metrical information about the scale of the simu- 
lation. Each subvolume is assumed to be homo- 
geneous and is provided with information con- 
cerning the numbers of simulation particles it 50 
contains (apportioned appropriately among 
chemical species and energy as required by the 
simulation parameters), its temperature (and 
pressure, if appropriate), its size and shape, and 
its location in the material energy system 60. It is 55 
also tagged with one or more material IDs de- 
scribing its composition and properties via a link 
to the material property data group 44 as descri- 
bed further below. The overall structure of the 



subvolumes can be static (determined at the be- 
ginning of the simulation process and then left 
unchanged) or dynamic (certain regions of sub- 
volumes are allowed to vary with time). An exam- 
ple of a simulation where dynamic subvolumes 
would be desirable is in simulation of chemical va- 
por deposition of thin films. The subvolumes can 
be regular or irregular in size, with the only con- 
straint that there can be no boundary mismatch 
between two adjacent subvolumes. The subvo- 
lumes can be monolithic, or separated into dis- 
tinct areas for convenience. For example, two 
subvolumes moving relative to each other can be 
used to simulate the effects of tribological inter- 
actions. As anot her example, in the laser heating 
embodiment (Figs. 9 and 10) the subvolumes can 
be conveniently separated into layers of different 
materials. Further, one or more individual layers 
can be configured to be patterned or composite 
layers by assigning different materials ID'S to in- 
dividual subvolumes in those layers. 

2. The materials property data group 44 contains 
basic physical and chemical information describ- 
ing the materials and chemical species present in 
the system 60. These can be organized as con- 
venient, with each grouping identified by an ap- 
propriate ID tag used for linking to the subvo- 
lumes and the event probability data group 46. In- 
formation which may be specified includes molar 
concentrations, heat capacity, thermal conduc- 
tivities, optical properties, densities, susceptibili- 
ties, emissivities, compressibilities, permeabili- 
ties and so on. Geometrical information about en- 
ergy and matter sources such as a laser beam or 
a gas jet would also be included. Each may have 
appropriate dependences on temperature, wave- 
length of light, time, pressure, chemical concen- 
trations, phase of matter, and local electric, mag- 
netic and gravitational fields. When new probabil- 
ities are to be updated, the event probability data 
group 46 is available from the data storage appa- 
ratus. A particular feature of this type of data 
structure is that the computational overhead of 
including very complex dependences is only 
slightly larger than including very simple depend- 
ences, or assuming that physical properties are 
constants. This allows realistic tracing of individ- 
ual events in different subvolumes to be made. 
The material property data group 44 can be inde- 
pendently generated and stored in the data stor- 
age apparatus and available for retrieval during 
the simulation and/or on-line control processes. 

3. The event probability data group 46 comprises 
subgroups of event probabilities organized ac- 
cording to a predefined order, and includes all 
events selected for simulation. It is linked to the 
subvolumes by assigning each subgroup mem- 
ber to a subvolume or a set of subvolumes. Prob- 
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abilities for each type of process can be listed 
separately in one or more subgroups. For exam- 
ple, in the laser heating embodiment the probabil- 
ities for energy absorption and heat flow are listed 5 
in separate arrays. Each subvolume is linked to 
one or more event probability data group mem- 
bers via its INDX. For example, in the laser heat- 
ing embodiment (Fig. 10), subvolumes and prob- 
abilities are organized into separate arrays, and 10 
INDX is a subvolume-array index. A specific 
probability array element can be accessed by us- 
ing an arithmetic expression to calculate the 
probability array index for that element from 
IMDX. In this manner, a simple connection be- is 
tween probability group member and subvolume 
or a set of subvolumes undergoing a particular 
process may be maintained. In other situations 
other methods of organizing the probabilities and 
defining INDX may be more useful. The advan- 20 
tage of providing separate event probability sub- 
groups for each distinct physical process is that 
it allows the explicit selection of the events to be 
simulated at the time the data structure 40 is con- 
structed. In this manner a very flexible, general- 25 
purpose simulator can be implemented which can 
be used to produce 3-D simulations of one or 
many simultaneous events. Because the prob- 
ability data is linked to the subvolumes 62 identi- 
fied by the associated subvolume INDXs, and 30 
each subvolume 62 is labelled with one or more 
material property IDs, the probability data group 
48 is automatically linked to "libraries" of informa- 
tion, concerning event processes and properties 
of the subvolumes defined by the material prop- 35 
erty data group 42 which are thus available for 
probability determination, 
i 4. The event process data group 48 contains all 

process functions, which may be a set of mathe^ 
matical representations or database tables or 40 
other kinds of relational definitions for process 
probability determinations. 
These process functions may be fixed, e.g. diffu- 
sion, or stored separately in a database before the 
simulation process, e.g. chemical reaction steps and 45 
kinetics. Typically these process functions will use va- 
rious material-specific parameters from the materials 
properties data group 44 and the system conditions 
as provided by the system state data group 42. The 
event process data group may be constructed to sim- 50 
ulate one or more of the following processes which 
may occur in the material-energy system 60: 

a. chemical transformation 

b. heat transfer 

c. mass transfer 55 

O diffusion 

O flow 

O convection 

O movement under an electric field 
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O movement under a magnetic field 
O movement under a gravitational field 
O movement under a centrifugal field 

d. absorption, reflection and emission of radiation 

e. compression by a static or time-varying force 
(stress or shock) 

f. phase transformation 

g. radioactive decay, transmutation 

h. mechanical and frictional forces 

The data structure 40 as shown in Fig. 3B there- 
fore has unique functional characteristics which are 
particularly useful for a spatially inhomogeneous sto- 
chastic simulator. Among these data groups, a linking 
relationship is provided by two common indexing 
parameters, i.e. the subvolume INDX and the mate- 
rial property ID. By the use of these common indexing 
parameters, a simulator is provided with information 
to first determine the types of physical or chemical 
events that will occur in different subvolumes 62 us- 
ing event probability data group 46. In response to the 
occurrence of an event in a subvolume 62 labelled 
with one or more material IDs, the simulator is also 
provided with material property data group 44 to up- 
date the system condition data group 42 and event 
probability data group 46 by use of the process func- 
tions contained in the event process data group 48. 

The data structure 40 thus provides to the simu- 
lator sufficient information to trace a three- 
dimensional time-progressive material energy sys- 
tem without explicitly and simultaneously managing 
the interactions between the subvolumes. The com- 
plexity of a spatially inhomogeneous stochastic sim- 
ulator is greatly reduced. Specifically, in randomly se- 
lecting and tracing an event in each time step for prop- 
agating the system, the simulator in the mean time 
also randomly selects a set of subvolumes 62 implic- 
itly. The effects caused by the spatially inhomogene- 
ous material distributions and system conditions 
among different subvolumes 62 are accounted for via 
the link INDX between the event probability data 
group 46 and the system state data group 42 and the 
link material ID between the system state data group 
42 and the material property data group 44. 

In addition to the advantage of reducing the spa- 
tial complexity of stochastic simulation, the data 
structure 40 also provides greater flexibility for the 
construction of a stochastic simulator. A stochastic 
simulator generally has a set of core operations which 
basically includes the steps of 1). random selection of 
a simulated triggering event and 2). propagating the 
system with the selected event based on a time step 
depending on the total probability. With the flexibility 
provided by the data structure40, this set of core sim- 
ulation operations once developed can be applied to 
simulate a wide variety of systems in a modularized 
fashion. 

Simulation of material-energy systems repre- 
sented by different sets of data groups can be conve- 
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niently accomplished by applying this core simulation 
module to different data structures 40 comprising 
these different sets of data groups. A more generic 
simulator is therefore constructed by coupling this 
core simulation module with the data structure 40 
thus enabling the simulator to be easily convertible 
from one type of software and hardware implementa- 
tion to another and/or for simulating different materi- 
al-energy systems. 

Fig. 5 shows the actual implementation of the 
data structure 40 in a software/hardware system 70. 
The data structure 40 including the system state data 
group 42, the material property data group 44, the 
event probability data group 46 and the event process 
data group 48 are stored in a data storage apparatus 
71. A processor system 72 including a simulating 
means 74 for simulating a material-energy system 60 
is used to perform the steps described in stage 2 and 
stages 3 of Appendix A which comprise: 

1. determining the initial probabilities of all the 
events defined in the event probability data group 
46; 

2. determining a time step as a function of the to- 
tal probability; 

3. normalizing the probability; such that the sum 
of all probabilities is equal to a predefined value; 

4. randomly selecting a simulation triggering 
event; 

5. with the occurrence of the selected triggering 
event, updating the system state data group 42 
which may comprise data of material concentra- 
tions, temperatures, pressures, etc, using the 
event process data group 48 and the data provid- 
ed by the material property group 44, and the ex- 
isting system condition provided by the system 
state data group 42; 

6. updating the event probability data group 46 by 
use of the material property data group 44 and 
the updated system state data group 42; and 

7. repeating the above steps 2-6 until the simu- 
lation process is completed. 

The simulating means 74 implemented on the 
processor system 72 can either be software or hard- 
ware or any combinatorial implementations. The proc- 
essor system 72 has a data bus 73 linked to the data 
storage apparatus 71 for receiving and sending data 
to the data storage apparatus 71 . The software/hard- 
ware system 70 further has a display apparatus 76 for 
display, and an user interface system 77 for sending 
commands to the software/hardware system 70 and 
a controlling means 78 for controlling the material-en- 
ergy system 60 which is connected through a 
data/command connecting line 79 to the controlling 
means 78. The controlling means 78 further can re- 
ceive sensor and measurement signals from the sys- 
tem 60 for updating the system state data group 42 
stored in the data storage apparatus 71 . Through the 
user interface system 77, an operator of the soft- 



ware/hardware system 70 can provide controlling 
parameters for the simulating means 74 such as the 
number of simulation steps, the numbers of particles 

5 to represent chemical concentrations and energy and 
the display interval for the display data which may in- 
clude spatial distributions of concentrations, temper- 
ature or pressure or other system parameters to be 
sent to the display apparatus. 

10 On-line real time control of the operation of sys- 

tem 60 is achievable by the software/hardware sys- 
tem 70 if the simulating means 74 can propagate the 
system performance of system 60 faster than the re- 
alistic physical propagations of the process events in 

15 the system 60. 

In other words, if the simulating means 74 can 
complete the simulation and prediction steps within a 
time span shorter than the time step used for system 
propagation, then a system operator or the controlling 

20 means 78 can use the system performance predic- 
tion to control the real time operation of system 60. 

The simulation process is further simplified in 
this invention to achieve higher simulation speed in 
order to achieve the purpose of a real time on-line 

25 control of the material energy system 60 with the soft- 
ware/hardware system 70. Higher speed simulation is 
achieved by reducing the stochastic selection proc- 
ess to a one-dimensional selection process by taking 
advantage of the relational-definition provided by the 

30 data structure 40 and by arranging the event prob- 
ability data group 46 in a predefined order. 

Instead of randomly selecting a triggering event 
in each subvolume and tracing these randomly se- 
lected events simultaneously in a spatially inhomoge- 

35 neous system as proposed by Gillespie, the sub- 
groups of the event probability data group 46 includ- 
ing all the events in all subvolumes arranged in a pre- 
defined fixed order are treated as a single entity. 
Since the probabilities of all the events in all the sub- 

40 volumes are summed to obtain the time step, treat- 
ment as a single entity allows a single time base to be 
established for the entire system. The value of each 
event probability may vary in each time-step caused 
by the changes of the system conditions in each sub- 

45 volume, but the 'stacking order' of the event probabil- 
ities as defined by the event probability array 46 is 
maintained unchanged. Propagation of the system 
simulation can be more expeditiously performed by a 
simplified process wherein for each time step the sim- 

50 ulating means 74 randomly selects one triggering 
event among the event probability data group 46 
comprising all event probabilities. This is achievable 
for a spatially inhomogeneous system 60 comprising 
many subvolumes wherein each subvolume may 

55 have different event probabilities, material concentra- 
tions, reaction rates, temperatures, pressures, etc, 
because there is a relational-definition, i.e. subvo- 
lume INDX, which provides a link between each event 
probability to one or more subvolumes 62 and thence 
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to a set of system characteristics as provided by the 
system state data group 42 and a set of material prop- 
erties through the material property IDs. 

Higher speed of system propagation is also ach- 5 
ieved because unlike the traditional stochastic simu- 
lators, the simulation means 74 in randomly selecting 
an event from the event probability data group 46 
does not process the system propagation based on a 
general assumption that the selected event would 10 
cause the conditions with their associated event 
probabilities in all the subvolumes to change. Instead, 
the selected event is associated with a particular set 
of subvolumes identified by a set of subvolume INDXs 
and the simulating means 74 has sufficient informa- 15 
tion to determine rf any additional subvolumes are af- 
fected by the selected triggering event. Updating of 
the system state data group 42 is only carried out for 
those limited number of subvolumes instead of the to- 
tal system as is generally performed in the conven- 20 
tional simulator. 

One specific example that demonstrates these 
advantages is illustrated in Figs. 6A, 6B and 7, where 
a system 80 (Fig. 7) simultaneously undergoing 
chemical reaction, mass transfer and heat flow. For 25 
simplicity of illustration, the system 80 undergoes the 
same three chemical reactions as that shown in Fig. 
2 and also incorporates mass and heat transfer within 
a 3x3x3 array of subvolumes (Fig. 7). The event prob- 
abilities for all possible processes in the system 80 30 
are collected in the set of sixteen probability arrays 
shown in Fig. 6A. These arrays comprise: 

1. the relative probabilities of each of the three 
chemical reactions occurring in a specific subvo- 
lume (3 reactions x 27 subvolumes for 81 terms), 35 

2. the relative probabilities of mass transfer of any 
of the four chemical species (A..D) in the x, y or 
z directions across any of the 18 suvolume inter- 
faces perpendicular to the direction (4 x 3 x 1 8 for 

216 terms), and 40 

3. the relative probabilities of heat transfer in the 
x, y or z directions across any of the 1 8 interfaces 
perpendicular to that direction (3 x 18 for 54 
terms). 

The categorization of processes into independent 45 
probability arrays in this manner offers great flexibil- 
ity in choosing which processes to include in a partic- 
ular simulation. There is a direct correspondence 
which uses the subvolume INDX between a particular 
subvolume, or interface between a pair of subvo- 50 
lumes, and a specific member of each of the sixteen 
probability arrays. In Fig. 7 this correspondence is il- 
lustrated for the particular case of transport in the z 
direction across an interface between two subvo- 
lumes. For example, the value of z-transport probabil- 55 
ity array member P[9] would represent the probability 
for transport from the subvolume at position x=2, y=2, 
2-1 to the subvolume at position x=2, y=2, z=2. Tak- 
ing (221) to be INDX for the source subvolume, and 



assuming P[9] is the probability of a selected event, 
INDX is determined from the number 9 using the cor- 
respondence equations given in the Figure. This par- 
ticular correspondence would apply to the arrays 
TAZ, TBZ, TCZ, TD2 and HZ in Figure 6A. Other cor- 
respondences would be established for other process 
probability data group members since different sub- 
volumes or subvolume interfaces would be involved. 

In the event selection process, rather than con- 
sidering each probability array separately as would 
be done following Bunker et a) and Gillespie, all six- 
teen probability arrays are collected sequentially and 
the union of these arrays is considered as a whole. 
After calculating the total probability of these sixteen 
arrays, a time step is calculated from the total prob- 
ability. A random number is generated and probability 
array elements are summed sequentially until a sub- 
total equal to or greater than the random number is 
found. The probability array element thus selected 
determines the event selected for the time step. This 
might be, for example, a member of the array TBY, 
corresponding to transfer of species B in the y-direc- 
tion between two subvolumes. Depending on the val- 
ue of the random number, any of the 351 elements of 
the sixteen probability arrays may be selected. The 
likelihood for selection of a given event is equal to its 
contribution to the total probability P. This is repre- 
sented by the lengths of the line segments in Fig. 6B. 

The flexibility and the modular organization pro- 
vided by the simulating means 74 and the data struc- 
ture 40 enable the software/hardware system 70 to 
be implemented in a wide variety of processors and 
data storage configurations. Taking advantage of the 
more expeditious simulation propagation process and 
the high speed processor systems with increased 
data storage capacity, the real-time control of an on- 
line system 60 as that shown in Fig. 5 is achievable. 

A preferred embodiment is described below as an 
example of such implementation. 

Fig. 8A shows a chemical reactor 300 comprising 
an inlet valve 302 for filling the reactor with a plurality 
of reagents, an outlet valve 304 for draining the reac- 
tor 300 after the consumption of the reagents reaches 
a predefined percentage, e.g. ninety-two percent, 
and a plurality of heaters 306 disposed on the bottom 
of the chemical reactor for heating the reagents con- 
tained therein. The reaction rates and the consump- 
tion of the reagents are dependent on the tempera- 
ture distribution in the reactor 300. The optimal level 
of reagent consumption, i.e. ninety-two percent, is to 
prevent excessive by-products from being produced. 
There are six temperature sensors 308 monitoring the 
temperature distribution of six subvolumes 31 0 of the 
reactor 300 and the temperature measurements from 
these six sensors 308 are converted as digital signals 
by an analog to digital converter 312 which are then 
received by an on-line process control simulator 320 
as the temperature input. The simulator 320 receives 
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the temperature signals as the initial condition before 
the simulation and periodically during the on-line con- 
trol process to periodically update a system state data 
group 322. The simulator also performs the simula- 
tion steps as described in Appendix A to predict the 
percentage of reagent consumption as a function of 
time. The simulator also predicts the time when the 
predefined optimal level of consumption is complet- 
ed. An output signal is generated before the optimal 
level consumption is reached to a chemical reactor 
controller 330 which automatically opens the outlet 
valve 304 at the predicted time to drain the products 
generated in the chemical reactor 300 and to initiate 
a refilling process to start the chemical reactor 300 
again. The use of a more complicated design of tem- 
perature control may be avoided while an optimal re- 
action completion level is maintained by use of the on- 
line simulator 320. A more flexible operation is also al- 
lowed because the chemical reactor can accommo- 
date a wider range of initial temperatures of the re- 
agents when the temperature distribution is moni- 
tored and controlled by the simulator 320 and the con- 
troller 330. 

Fig. 8B shows another chemical reactor system 
400 comprising a chemical reactor, a simulator and a 
reactor controller similar to the chemical reactor sys- 
tem 300 as shown in Fig. 8A. The chemical reactor 
system 400 differs from the reactor system 300 that 
in the reactor system 400 the heaters 408 are control- 
led by the simulator 420 and the reactor controller 430 
to achieve a certain reaction rate. After receiving the 
temperature signals from the temperature sensors 
412, the simulator 420 not only updates the system 
state data group, it also predicts a temperature profile 
to achieve a specific reaction rate based on the ex- 
isting temperature distribution. The predicted temper- 
ature profile is transmitted to the reactor controller to 
adjust the heaters 408 whereby the reaction rate and 
the production time of the chemical reactor 400 can 
be actively controlled. By adding the heater control 
mechanism to the reactor controller, a more precise 
control of the production time is achieved. 

Fig. 9 shows a material-energy system 500 
wherein a laser beam 502 is irradiated on the top sur- 
face of the system 500 which generates an energy in- 
put, i.e. rate of heat, Q, flow into the system. The heat 
flow Q is a function of the temperature-dependent 
optical properties of the cell onto which the laser 
beam 502 is incident. The power of the incident laser 
beam 502 is represented by I and the temperature- 
dependent reflectivity is represented by R{T) where T 
is the temperature. The rate of heat transfer to an ad- 
jacent cell depends on the temperature-dependent 
thermal conductivity K(T), the temperature difference 
between the first and second subvolumes, the con- 
tact area A between the subvolumes and the distance 
between their centers, d. The instantaneous heat 
content Q of a subvolume determines its temperature 
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according to the relation T=Q/C[p]M, where C[p] is 
the heat capacity and M is the number of moles of ma- 
terial in the subvolume. 

5 Fig. 10 shows the material-energy system 500 

which is divided into a plurality of subvolumes 510. 
The subvolumes 510 are cyiindrically symmetric to al- 
low two-dimensional coordinates to be used. The sys- 
tem 500 comprises three different layers of materi- 

10 als,i.e. layer 520, layer 522 and layer 524 which may 
in turn comprise one or more slices where the thick- 
ness of a slice is the thickness of the subvolume. All 
light is absorbed in the top slice 512. The rate of en- 
ergy variations in the subvolumes near center of the 

15 top slice 512 is expected to be greater than the sub- 
volumes 510 which are remote from the laser beam 
502 which is the heat source. The material property 
data group provides information for R(T) which is the 
temperature-dependent reflectivity of the top slice 

20 512 while the thermal conductivity K(T), density 
&rho., and heat capacity C[p] are provided for each 
of the three layers 520, 522 and 524. The simulator 
predicts the temperature distribution as a function of 
time. The material-energy system 500 can be a thin 

25 film comprising three layers of materials deposited on 
a substrate (not shown) whereupon a laser beam 502 
is applied to the top slice 512. By the use of a simu- 
lator 540, a laser source 550 which generates the las- 
er beam 502 can be controlled to produce a prede- 

30 fined temperature profile in the three layers 520, 522 
and 524 of in the thin-film system 500. 

By examining the simulation results as that 
shown in Figs. 11,12 and 13, the temperature profile 
in the system 500 can be accurately controlled ac- 

35 cording to the preferred embodiment as that shown 
in Fig. 10. Fig. 11 is a comparison of the results of a 
temperature calculation performed with simulator 
540 to the results obtained by an analytic solution to 
the heat-diffusion equation for pulsed infrared laser 

40 heating of silver. The analytical solution is disclosed 
in a publication by J. M. Hicks, etal., Physical Review 
Letters, Vol. 61, No. 22, pp. 2588-2591, 1988. The 
materials and beam parameters used in simulator 540 
are the same as those used in the analytical solution. 

45 The results plotted are the time-dependent tempera- 
tures of the center subvolume of the top slice 512, 
onto which the laser beam 502 is entirely incident and 
absorbed. Its radius is 0.9 mm and its thickness is 27 
nm. 

50 Figure 12 is a comparison of results of tempera- 

ture as predicted by simulator 540 to an analytic sol- 
ution for cw visible laser heating of silicon. The ana- 
lytical solution is disclosed in a publication by E. Lar- 
iokapis et al., Journal of Applied Physics, Vol. 57, No. 

55 12, pp. 5123-4126, 1985. Both solutions utilize the 
temperature dependent thermal and optical proper- 
ties of silicon, and the temperature distributions at the 
top slice 512, i.e. the silicon surface, are compared. 
Results for several spatial geometries for the stochas- 

9 
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tic simulation are shown. 

Fig. 13 shows the results of the temperature dis- 
tribution predictions for cw laser heating of 2 &mu.m 
thick, 25 - 50 &mu.m diameter films of Cr, Mo and W 5 
on large silicon wafers. Results for the top slice 512 
of each metal at the various laser powers (in MW/cm 2 ) 
are shown. The dashed vertical lines delineating re- 
gions of pure metal and metal contaminated by car- 
bon and oxygen are for comparison to experimental 10 
data, and serve to determine a minimum temperature 
for obtaining clean metal films in laser chemical vapor 
deposition from the Cr, Mo and W hexacarbonyls. The 
optical properties of the metals were assumed to be 
temperature-independent, all thermal properties for 15 
the metals and the silicon substrates were tempera- 
ture dependent. 
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1 . A simulator system for simulating and/or control- 
ling a time-progressive material-energy system 
comprising at least one data storage apparatus 
including (F1) a data correlating means; a data 25 
structure (40) stored in said data storage appara- 
tus divided into data groups defining geometric 
configuration, system state, material properties, 
material-energy processes, and probabilities of 
occurrence of said material-energy processes; 30 
said data correlating means being capable of de- 
fining at least one functionai-relationship link be- 
tween data stored in said data groups. 

2. A system as claimed in Claim 1 , further compris- 35 
ing at least one processor (72) including a data in- 
put means for receiving said data from said data 
storage apparatus, the processor further having 
means (74) for performing a simulating process 
utilizing said functional-relationship link for corre- 40 
lating among said data comprising the steps of: 

(a) computing an event probability for each of 
the possible said material-energy processes 
and determining a time step by using a total 
probability of occurrence by summing up all of 45 
said material-energy probabilities in said geo- 
metric configuration; 

(b) randomly selecting a material-energy 
process event for said time step and using 

said material property data and said geomet- 50 
ric configuration data for updating said system 
state data resulting from the occurrence of 
said selected material-energy process event; 

(c) repeating said steps (a) and (b) until a pre- 
defined simulation end time is reached; and 55 
the processor further having a control means 

(78) capable of utilizing said updated system 
state as obtained from step (b) for controlling 
said time-progressive material-energy sys- 

10 



tern. 



3. A system, as claimed in Claim 2 wherein the step 
of randomly progressively summing the probabil- 
ities of occurrence of all processes in a fixed or- 
der and comparing the sum to a random number. 

4. A system as claimed in any of claims 1 to 3 where- 
in said data storage apparatus further includes a 
relational database for storing said data and said 
data correlating means includes a relational da- 
tabase manager. 

5. A system as claimed in any preceding claim fur- 
ther comprising a display means for displaying 
said data whereby an operator of said material- 
energy system can visualize said displayed data. 

6. A system as claimed in Claim 2 or Claim 3 in 
which the system state data group includes a 
plurality of subvolumes each being identified by 
a subvolume index (INDX), the system state data 
group including a plurality of data subgroups for 
defining a system state for each of said subvo- 
lumes, the material properties data group in- 
cludes a plurality of data subgroups each defin- 
ing a plurality of materials and a plurality of ma- 
terial properties for each of said materials con- 
tained in said subvolumes identified by said sub- 
volume INDX, the material energy data group de- 
fines a plurality of events each having a function 
for defining a time dependence of said event as 
a function of said system state in each of said sub- 
volumes identified by said subvolume INDX; and 
the probabilities data group defines a probability 
of occurrence for each of said events in each of 
said subvolumes depending on said system state 
and said material property in said subvolume. 

7. The simulator system of claim 12 wherein said 
processor further comprises at least one data 
port for receiving system state data from said ma- 
terial-energy system for updating said system 
state data group in at least one of said time steps. 

8. A system as claimed in Claim 7 wherein the time- 
progressive material energy group defines a mul- 
ti-layer thin film formed on a top surface of a sub- 
strate with a laser source irradiating a laser beam 
onto said thin film; said data port of said proces- 
sor furt her includes at least one temperature sen- 
sor disposed in at least one of said subvolumes 
for measuring a temperature in said subvolume; 
and said control means of said processor further 
controls said laser source to generate a prede- 
fined temperature profile in said multi-layer thin 
film. 
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9. A system as claimed in Claim 7 wherein the ma- 
terial energy group defines a chemical reactor 
having a reactor chamber, an inlet valve, an outlet 
valve and a plurality of heaters in said reactor 5 
chamber wherein a plurality of reagents are re- 
ceived from said inlet valve into said reactor 
chamber for making at least one product to be 
output from said outlet valve at a predefined sys- 
tem state; said data port of said processor further 10 
includes at least one temperature sensor dis- 
posed in at least one of said subvolumes for 
measuring a temperature in said subvolume in 

said reactor chamber; and said control means of 
said processor further controls said heaters, said is 
inlet valve and said outlet valve to generate a pre- 
defined operation condition in said reactor cham- 
ber. 

1 0. A method of providing a data structure for a spa- 20 
tially inhomogeneous stochastic simulator com- 
prising the steps of generating system state data 
including a plurality of subvolumes each identi- 
fied by a subvolume index (INDX), subvolume 
system state data and at least one material prop- 25 
erty identification (ID); generating material prop- 
erty data including the properties of the materials 
contained in said plurality of subvolumes; gener- 
ating event process data for def ining a plurality of 
events each having a function for defining a time 30 
dependence of said event as a function of said 
system state for each said subvolume identified 

by said subvolume INDX; generating event prob- 
ability data for defining a probability of occur- 
rence for each of said events in each said subvo- 35 
lume depending on said system state and said 
material property in said subvolume; and storing 
said system state data, said material property 
data, said event process data and said event 
probability data in a data storage apparatus 40 
wherein said subvolume INDXs and said material 
property identification (ID) provide a correspon- 
dence between all of said data. 

11. A method as claimed in claim 10 for simulating a 45 
time-dependent material-energy system by using 

an executable program on a computer comprising 
the steps of: 

(a) generating a system state data group in- 
cluding a plurality of subvolumes each identi- so 
fied by a subvolume index (INDX), system 
state data subgroup and at least one material 
property identification (ID); 

(b) generating a material property data group 
including a plurality of data subgroups for de- 55 
fining the material properties contained in 
said subvolumes identified by said associated 
subvolume INDX; 

(c) generating an event process data group for 



defining a plurality of events each having a 
function for defining a time dependence of 
said event as a function of said system state 
for each said subvolume identified by said 
subvolume INDX; 

(d) generating an event probability data group 
for defining a probability of occurrence for 
each of said events in each of said subvo- 
lumes depending on said system state and 
said material property in said subvolumes; 

(e) storing said system state data group, said 
material property data group, said event proc- 
ess data group and said event probability data 
group in a data storage apparatus wherein 
said subvolume INDXs and said material 
property identification (ID) provide a corre- 
spondence between all of said data groups; 

(f) computing said event probability for each of 
said events in each of said subvolumes and 
determining a time step by using a total prob- 
ability of occurrence by summing up all of said 
event probabilities; 

(g) randomly selecting an event for said time 
step and using said material property data 
group for each of said associated subvolumes 
for updating said system state data group rep- 
resenting the system state resulting from the 
occurrence of said selected event in each of 
said subvolumes; and 

(h) repeating said steps (f) and (g) until a pre- 
defined simulation end time is reached. 
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Event Selection for Homogeneous Chemical Reaction 



FIG. 2A 
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Select Event in Simulation of Chemical Reaction, Mass 
Transfer, and Heat Row in Cube (3x3x3 cube, 3 reaction 
stepsj 
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7 ' (18 S demenf) r ' ^ B ' X ° ?redl ° n: ™11H.211- 
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Example: Probabilites for xyz 
transfer in z direction 



x'yY cell pairs, 



Pi 

P3 

P5 

P? 

P9 

P11 

p 13 

P15 

P17 



P2 

P4 

P6 

P8 

PlO 

P12 

PH. 

P16 

P18 



111~112~113 
121 «* 122«123 
131~132~133 
211~ 212^213 
221~222~223 
231~232~233 
311~312~313 
321~322~323 
331~332~333 



FOR PROBABILITY Pn TO TRANSFER FROM XYZ XY (2+1) 
(WITHIN UMITS OF VOLUME 80) 1 ' 

X= [(n-1) div 6] +1 

y= [[(n-1) div 2] mod 3] +1 

z= [(n-1) mod 2] +1 
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(54) Stochastic simulation system. 

(57) The present invention comprises a simulator 
including at least one data storage apparatus 
(71) and at least one processor (72). The 
simulator further has a material property data 
group (44), a system state data group (42) com- 
prising geometrical data for defining a plurality 
of geometrical subvolumes and system condi- 
tion data representing the state of each sub- 
volume of a simulated system, an event 
probability data group (46) including a plurality 
of event subgroups for uniquely identifying a 
set of events each having a probability of oc- 
currence represented by the event probability, 
each of the event probabilities being associated 
with at least one of the geometrical sub- 
volumes, and an event process data group (48) 
including the time dependent process functions 
for all events which can occur in the system. 
These data groups and at least one functional- 
relationship link for correlating all of them, 
whereby the material energy system can be 
represented by such data groups without re- 
quiring an explicit multiple dimension data 
structure, are stored in the data storage ap- 
paratus. The processor includes a data input 
means (73) to read these data groups from the 
data storage apparatus. The processor further 
has means for performing a simulating process 
comprising the steps of: (a) computing the 
probability of occurrence for each of the events 
and determining a time step by using a total 
probability of occurrence as a sum of the prob- 
ability of occurrence for each of the events ; (b) 
randomly selecting a simulated event for the 
time step and updating the system condition 



data group representing the system state resul- 
ting from the occurrence of the selected simu- 
lated event; and (c) repeating the steps (a) and 
(b) until a predefined simulation end time is 
reached. 
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